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We discuss a more general class of phantom (p < —g) cosmologies with various forms of both 
phantom {w < —1), and standard {w > —1) matter. We show that many types of evolution 
which include both Big-Bang and Big-Rip singularities are admitted and give explicit examples. 
Among some interesting models, there exist non-singular oscillating (or "bounce") cosmologies, 
which appear due to a competition between positive and negative pressure of variety of matter 
content. From the point of view of the current observations the most interesting cosmologies are 
the ones which start with a Big-Bang and terminate at a Big-Rip. A related consequence of having 
a possibility of two types of singularities is that there exists an unstable static universe approached 
by the two asymptotic models - one of them reaches Big-Bang, and another reaches Big-Rip. We 
also give explicit relations between density parameters fl and the dynamical characteristics for 
these generalized phantom models, including higher-order observational characteristics such as jerk 
and "kerk". Finally, we discuss the observational quantities such as luminosity distance, angular 
diameter, and source counts, both in series expansion and explicitly, for phantom models. Our 
series expansion formulas for the luminosity distance and the apparent magnitude go as far as to 
the fourth-order in redshift z term, which includes explicitly not only the jerk, but also the "kerk" 
(or "snap") which may serve as an indicator of the curvature of the universe. 

PACS numbers: 04.20. Jb,98. 80. Jk,98.80.Cq 



I. INTRODUCTION 

Phantom is matter which possesses a very strong nega- 
tive pressure and so it violates the null energy condition 
g -\- p > 0, where g is the energy density and p is the 
pressure. Consequently it violates all the remained en- 
ergy conditions: the strong, the weak and the dominant. 
In terms of the barotropic equation of state p = wg {w = 
const.) phantom is described as matter with w < —I. 
The story of phantom and its observational support is 
not quite analogous to the story of quintessence or dark 
energy - the matter which violates the strong energy con- 
dition g-\-3p > 0, i.e., the matter which allows w < —1/3. 
Namely, despite its theoretical investigations P, , and 
even admission for the early universe inflation, it was 
hardly accepted as reality before the first supernovae 
data was revealed A different story happens with 
phantom. A more detailed check of the data suggests 
that the current supernovae not only gives evidence for 
quintessence with —1 < w < —1/3, but also admits phan- 
tom with w < — 1 at a high confidence level 0|. The point 
is that in view of the data from supernovae there is no 
sharp cut-off of the admissible values of w at w = —1. 
Instead, the data even favours the values of w which are 
less than minus one. 

Phantom matter was first investigated in the current 
cosmological context by Caldwell 01 , who also suggested 
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the name referring to its curiosity. In fact, phantom (or 
ghost) must possess negative energy which leads to insta- 
bilities on both classical and quantum level Q ■ Since 
it violates all possible energy conditions, it also puts in 
doubt the pillars of relativity and cosmology such as: the 
positive mass theorems, the laws of black hole thermody- 
namics, the cosmic censorship, and causality in relativity 
theory H, In this context phantom becomes a real 
challenge for the theory, if its support from the the super- 
novae data is really so firm. The thing is that the main 
support for the phantom comes from supernovae dim- 
ming rather than from large-scale structure and cosmic 
microwave background observations. That is why some 
alternative explanation of the dimming of supernovae due 
to a conversion of cosmological term and photons into ax- 
ions was made 10 . From the theoretical point of view - a 
release of the assumption of an analytic equation of state 
which relates energy density and pressure and does not 
lead to energy conditions violation (except for the domi- 
nant one) may also be useful 11]. As for the explanation 
of the supernovae data phantom is also useful in killing 
the doubled positive pressure contribution in braneworld 
models 

Phantom type of matter was also implicitly su gges ted 
in cosmological models with a particle production [13 , in 
higher-order theories of gravity models , Brans-Dicke 
models for w < — | , in non- minimally coupled scalar field 
theori es [l5l , in mirage cosmology of the braneworld sce- 
nario |l6|. and in kinematically-drivcn quintessence (k- 
essence) models example. The interest in phan- 

tom has grown vastly during the last two years and var- 
ious aspects of phantom models have been investigated 
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One of the most interesting features of phantom models 
is that they allow for a Big- Rip (BR) curvature singular- 
ity, which appears as a result of having the infinite values 
of the scale factor a{t) ^ cx) at a finite future. However, 
as it was already mentioned, the evidence for phantom 
from observations is mainly based on the assumption of 
the barotropic equation of state which tightly constraints 
the energy density and the pressure. It is puzzling [Tll | 
that for Friedmann cosmological models which do not 
admit an equation of state which links the energy den- 
sity g and the pressure p, a sudden future singularity of 
pressure may appear. This is a singularity of pressure 
only, with finite energy density which has an interesting 
analogy with singularities which appear in some inhomo- 
geneous models of the universe Q| ■ 

Recently, phantom cosmologies which lead to a 
quadratic polynomial in canonical Friedmann equation 
have been investigated However, only the two 

types of phantom matter with the equations of state 
p = — (4/3)f) (phantom) and p — —(5/3)^3 (superphan- 
tom) were admitted. In this paper we extend our dis- 
cussion into four phantom fluids by adding subphantom 
(p — —2g), and hyperphantom {p — —(7/3)g), and con- 
sider all types of matter with —7/3 < w < 1 which 
completes a general discussion of the Friedmann equa- 
tion with many-fluid matter source given in series of pa- 
pers ;55]. As it was found in there exist interesting 
dualities between phantom and ordinary matter models 
which are similar to dualities in superstring cosmologies 
These dualities were generalized to non-flat and 
scalar field models ^8t> i^j i brane mod els i4Ll. and found 
to also be related to ekpyrotic models [S^ ISJ . 

In Section II we write down the field equations for 
Friedmann cosmology which involve all standard and 
phantom matter fluids and express them in terms of the 
observational parameters such as the dimensionless en- 
ergy density parameters fl, Hubble parameter H, decel- 
eration parameter q, and the higher derivative param- 
eters like jerk j, "kerk" k etc. [Hi El |H| which may 
serve as indicators of the equation of state (statefinders) 
and the curvature of the universe. We also give explicit 
expressions relating these quantities between themselves, 
which allows to eliminate some of them, when the com- 
parison with astronomical data is required. Some inter- 
esting symmetries in these relations, which refer to the 
appearance of the density parameters, are very clearly 
visible. 

In Section HI we study some exact cosmological mod- 
els which lead to a polynomial of the third and the fourth 
order in the canonical Friedmann equation. This is the 
next step to what has been done in our earlier reference 
[isf. since now the solutions are given in terms of the 
elliptic functions only. Because of the tenth degree poly- 
nomial in the basic canonical equation which involves all 
the fluids from w = —7/3 to w = 1, we discuss only the 
models which, after some change of variables, can be re- 



duced to the case of the third or fourth degree canonical 
equation polynomial. All the possible cosmological mod- 
els which are allowed in such a general case are discussed 
in the two appendices A and B. 

In Section IV we present the discussion of the obser- 
vational quantities (luminosity distance, angular diame- 
ter, source counts) both in terms of the series expansion 
which involves all the fluids under studies together with 
the higher-derivative parameters, and in some cases ex- 
plicitly given by elliptic functions. Similar results were 
obtained for quintessence (— 1 < ui < 0) models in an ar- 
chaic reference . In Section V we give the conclusions. 



II. FRIEDMANN EQUATION WITH 
PHANTOM AND OBSERVABLE PARAMETERS 

Phantom is a new type of cosmological fluid of a very 

strong negative pressure which violates the null energy 
condition (NEC) Hi] 

p + g>0, (H.l) 

i.e., it obeys a barotropic equation of state 

P=h-'^)Q = wg, {11.2) 

with negative barotropic index 

j^w + KO , (H.3) 

or, as more commonly expressed nowadays, with 

w < -1 . (H.4) 

Phantom also violates the strong, weak, and dominant 
energy conditions (SEC, WEC, and DEC, respectively) 

g + p > 0, and g + 3p > , (II. 5) 

g > 0, and g + p>0, (II. 6) 

e > 0, and - g < p < g . (II. 7) 

Phantom may be included into the basic system of 
equations for isotropic and homogeneous Friedmann uni- 
verse which reads as 

K b? 

K^g = -A + 3^ + 3^, (H.8) 

K^p = A- 2 (11.9) 

a a 

where a{t) is the scale factor, K ^ ±1 the curvature 
index, A the cosmological constant, n"^ ~ 8ttG - the Ein- 
stein constant. 

Using equations (jll.8|) and (|II.9|I . one gets the Fried- 
mann equation in the form 

a' ^ A 



3 



After imposition of the conservation law 
ga^'t = (3/k2)Ct, = const. , 



one gets Eq. (jII.10|) as follows 

2 
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(11.11) 



(11.12) 



The standard types of cosmological matter are: stiff-fluid 
(7 = 2), cosmic gas (7 = 5/3), radiation (7 = 4/3), 
dust (7 = 1), cosmic strings (7 = 2/3), domain walls 
(7 = 1/3), and cosmological constant (7 = 0) and the 
models which involve them have been studied in detail |2, 
Issf . They, in fact, admit a large class of oscillating (non- 
singular) solutions due to a balance between the positive 
and negative pressure domination periods. This is in full 
analogy to what happens in ekpyrotic scenario [sOj , where 
the positive potential energy dominates throughout some 
period of the evolution (including now) , and the negative 
potential energy dominates near the bounce. 

In this paper, which follows our earlier Ref. |45l |. we 
extend the discussion of the models onto the case of neg- 
ative barotropic index phantom matter with 7 = —1/3 
(phantom), 7 = —2/3 (superphantom) , 7 = — 1 (sub- 
phantom), and 7 = —4 /3 (hyp erphantom) . Due to a 
phantom duality psl liM lig . l51| 
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these phantom types of matter are dual to domain walls 
(7 = 1/3), cosmic strings (7 = 2/3), dust (7 ~ 1) and 
radiation (7 — 4/3), respectively. 

The simplest way to consider these dualities is to in- 
tegrate the Friedmann equation |^ 
fluid only, which gives 



for one general 



a(t) = ao I i 
for ordinary (7 > 0) matter, and 

a{t) ^ao\t 1"^ 



(11.14) 



(11.15) 



for phantom 7 = — | 7 |< 0. In both cases there is a 
curvature singularity at i = 0, but in the former case it 
is of a Big-Bang type, while in the latter case it is of a 
Big-Rip type. From the current observations, however, 
it is reasonable to choose the solution (|II.14I) for positive 
times t > 0, and the solution H11.15() for negative times 
t < 0. In this context the age of the universe in a stan- 
dard matter case is the time from Big-Bang at t = 0, to 
the present time t = to, and reads as 
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(11.16) 



while in a phantom case one can only calculate the time 
from the present t = —to to a Big-Rip at t = which 
reads as 



to 
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Similar investigations should be taken into account, if one 
wants to calculate the redshift of an observed galaxy. If 
a light ray was emitted at the time ti in the past, then 
in the standard matter case, the redshift reads as 



1 



a{ti) 



(11.18) 



We can see that the redshift tends to infinity, if ti ap- 
proaches zero (i.e. at Big-Bang). The situation is differ- 
ent for phantom, since the redshift is given by 



ajto) 
ait,) 



(11.19) 



and the redshift tends to infinity if the time ti approaches 
minus infinity. On the other hand, if the time of emission 
of the light ray is the time of Big-Rip, which is at ti — 0, 
then we have formally z = —1. 

Let us now define the appropriate cosmological param- 
eters which characterize the dynamics of the evolution of 
the universe which are: 
the Hubble parameter 



a 

the deceleration parameter 

1 a 
' = -lPa = 
the jerk parameter [s^ 



aa 
'7^ 



J = 



1 d 



aa 



(11.20) 



(11.21) 



(11.22) 



and the "kerk" (snap j53!|) parameter (which includes the 
fourth derivative of the scale factor a - we can carry on 
to call the higher derivative parameters "lerk", "merk", 
"nerk", "oerk" etc., respectively) 



aa 



(11.23) 



a a- 

The application of the definitions of the parameters 
and ljll.23p gives the following 



H = -H^iq + 1) , (11.24) 

ii = H^{j + 3q + 2) , (11.25) 

H = -H*[k + 4:j + 3q{q + A) + 6] . (11.26) 

A comparison of phantom models with observational 
data requires the introduction of dimensionless density 
parameters 



n 



KO 



1^A„ = 



3Hi 
K 

Aq 



(11.27) 
(11.28) 
(11.29) 
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where x = hp, bp, sp, ph, w, s,m, r, g, st (hyperphantom, 
subphantom, superphantom, phantom, domain walls, 
cosmic strings, dust, radiation, cosmic gas, and stiff-fluid, 
respectively), and the index "0" means that we take the 
parameters at the present moment of the evolution t — Iq. 

Applying the relation (K' = Cs — K, where Cg — 
for 7 = 2/3, cf. CniJ) 



Q.K'0 — ^ 



KO, 



(11.30) 



we realize that for a = aq, the Friedmann equation ljII.12|) 
can be written down in the form 



not appear in ljll.36|l . Next step is to express the kerk 



parameter by 



= qojo + GOnsto + SOflgo + U^ro + 3f} 



ifipO 



(11.37) 



Here, in turn, none of the three parameters: cosmic string 
energy density parameter Qgo, domain wall energy den- 
sity parameter Qwo, and the cosmological term VIa^ does 
not appear explicitly in ljll.37p . although they are in- 
cluded in both qq and jo- 



^hpQ + ^bpO + ^spO + ^phO + ^Ao (11.31) 



From the field equations ljll.8|) and ljll.9|) we can also 
derive the following relations between the observational 
parameters 



K 
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(11.32) 



spO 



^bpQ — 2flhpO — qo — 1 



K' 5 3 



— f2u,o ~ 2^p'* 



(11.33) 



^spo — 2^fcpo ~ 2flhpa — (?o ~ 1 , 



and 



— ^ — 2flstO + -r^^gO + ^rO + -^^'mO ^ QO ^ -^^wO 



3 5 

-^^phO ^ 2QspO ^ -^^bpO ^ SflfipO ■ (11.34) 



3H, 



Alternatively, one can use pi.32|l and ljll.34(l to express 
the deceleration parameter go as 

3 11 
go = 2risto + 2^90 + ^ro + 2^»"0 ^ 2^™" (11.35) 

3 5 
~ ^Ao — 2^p'»o ~ 2r2<,po ~ 2^''po ~ 3f2/ipo • 

Note that cosmic string energy density parameter flso 
does not appear in pi. 3 5(1 . Similarly, it is possible to 
express the jerk parameter 

jo = lOr^stO + 6r2gO + 3QrO + ^mO + ^Ag 

+ Sflpho + Q^spo + lOilbpo + 15r2;ipo ■ (n.36) 



III. EXAMPLES OF PHANTOM FRIEDMANN 
COSMOLOGIES 



Using the new variables 



, T = Hot, 



(IILl) 



one turns the basic Friedmann equation 1)11. 12p into the 
form 



^hpoy^ + ^bpoy^ + ^spoy* 



+^phoy^ + ^Aoy'^ + ^woy + ^K'o (ni.2) 
+^moy~^ + ^roy^^ + ^goy^^ + ^stoy^'^- 

The discussion now refers to an investigation of the 
cosmological models which are implied by this Friedmann 
Of course, the integration of a general 



equation (|III.2|) 
case would require the application of abelian functions, 
since we have the tenth degree polynomial to the right- 
hand side of 

-I.2II can be reduced to the 



(|III.2|) . In this paper we will discuss only 
the cases in which Eq. 
equation with the third or fourth degree polynomial in 
the canonical equation. Phantom cosmologies which re- 
sult from the second degree polynomial were discussed 
in |4^ . Since a general mathematical discussion of the 
canonical equations is quite extensive, we put most of it 
into the Appendices A and B, leaving only some examples 
of solutions in the main body of the paper. One example 
is out of the category of the third degree polynomial in 
the canonical equation, and another one is out of the cat- 
egory of the fourth degree polynomial. The last example 
is out of the category of the second degree polynomial 

(cf. H). 



A. Superphantom, A-term, strings, and radiation 
models 

The simplest case which can be integrable in terms 
of elliptic functions is when we neglect all the fluids ex- 
cept for superphantom, A-term, strings, and radiation in 



In fact, neither cosmic string energy density parameter 
risOj nor domain wall energy density parameter flwo, do 



^ ) = rispoy'^ + ^Aoy^ + ^k'o + ^roy ^ ■ (iii-3) 
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Since we deal with even powers of y only, we can intro- 
duce a new variable 



x = y , 

so that the equation (|111.3|l reads as 



dx\ 3 

— I = ilspQX 



^KOX + ^K'QX + il 



rO 



(III.4) 



W{x) , 



(ni.5) 

and this is an example of the third degree polynomial 
W{x) in canonical equation discussed in the Appendix 
A. From the definition 1)11.27(1 we conclude that ilspo is 
always greater than zero, and so we deal with the ap- 
propriate solutions of section A.l. of Appendix A. The 
discriminant A of W{x) is, in turn, a polynomial of the 
second order in flgpo with a negative coefficient in front 
of ^Ipo- Its roots are 



(III.6) 

The sign of A depends on the sign of r^'^/g — Sr^T-o^Ao- 
For A, in order to be positive, we must have 

^/c'o 3r2,oriAOj r^sp- < ^spo < ^sp+ ■ 

(III.7) 

The existence of the appropriate values of the parameters 
follows from a general analysis of the behaviour of the 
polynomial W{x). For A > 0, the solutions are given by 
(fO)l and l|X7|l . i.e.. 



y- 
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p(2t - To -f W3) 



AO 



spa 



P(2t - To) 



AO 



317 



spQ 



with the invariants 



92 



AO 



-ilgpO^K'O, 



12 

1 1 3 1 2 

53 = -^^spQ^AO^K'O - ^^AO ^ J^^spQ^rO- 

The solution y_ is oscillating and non-singular (ONS) 
type, while the solution is singular and evolves from 
a Big- Rip to a Big- Rip (BR-BR type) (compare Fig. 

Because of the change of the variable, 2r appears as 
an argument, and the "physical" period is not 2wi, but 
(jji. For the solution y- to be possible, the roots of W 
must be such that 63 < 1 < 62. This condition can be 
simplified by analyzing the position of the minimum of 
W. It yields 

1 3 

When the opposite is true, the other solution is valid. 
Moreover, as W{0) = Vlra > 0, we have 63 < 0, so that 
the oscillating solution y+ necessarily possesses the y = 



singularity. For the j/_ solution, we obtain the singular 
behaviour when QspO^K'o > and > 0. Otherwise 
ei > 0. 

If A is to be equal to zero (case A.l.b. of Appendix 
A), its discriminant must be non- negative, i.e., 

^K'o ^ Silrof^Ao, and flspo = ilsp± ■ (III.8) 
On the other hand, the root of W{x) is triple, if 



A — i7^Q — 3r2spo^A''o — 



(III.9) 



and 



so that the solutions are 



~2 

Vo 



and 



rO 



\ ^spO 



VlrO 



^spO f7spo(2T — To) 



(III.IO) 



Here the tilde indicates, that the Hubble constant Hq = 
0, so that the normalization is not possible. As the root 
is non-positive (VF(0) > 0), the solution yo is both BB 
and BR singular, whereas the solution yo is not physical 
(compare Fig. 0J. 

Another possibility is a double root. Both of the con- 
ditions piLlOjl . must be broken then (if only one was, 
the discriminant would not be equal to zero). There are 
two subcases, depending on the sign of 173. 

With 53 > 0, we have a stable static (SS) Einstein 
universe y = 0, if VIk'o = ^r o = 0, and another model of 
the universe given by 



Vl = ^ <! 2 



3 tan 



spO 



\/93{t - To) 



17 AO 



^^spO 

(IH.ll) 

which can be of BR-BR type, if f^Ao < 0, and BB-BR 
type, if r^Ao > 0. 

For 33 < the solutions are exactly the same as in 
section A.l.b. of Appendix A, i.e.. 




AO 



317spo 

17ao 
3Hsp0 

(III.12) 



The first solution is non-singular, if a static solution y — 
exists. The second is BB singular, and the third one 
might be free from the singularity, if the double root < 
Vo < 1- 

Finally, A is negative, when Vlspo < Vl^p- or Vlspo > 
Vlsp+. For Vlj^iQ < SflrQ^lAo, this is the case for all values 



6 



of i^spo- The solution is that of Appendix A. I.e. and, as 
above, W{0) > 0, which means it is not only BR singular, 
but also BB singular, despite the case presented in Fig. 
|S1 The exact solution is 



Vo 



-p(2t - To) - 



AO 



with "period" being 2TZe{uJi). 



(III.13) 



B. Hyperphantom, A-term, and radiation models 



The Firedmann equation (jIII.2|) now reduces to 



dy 
du 



(III. 14) 



We use the same variables as in the superphantom case 



x=y , 

so that (|III.14|I reduces to 

2 



(III.15) 



(III.16) 



which can be easily analyzed, as it is biquadratic. This 
is an example of the fourth degree polynomial W{x) in 
the canonical equation discussed in the Appendix B. The 
four roots are given by 



2n 



hpO 



(III.17) 



where JIao = 1 — ^hpo ~ ^ro-, and we can analyze the 
subcases with respect to the {i^hpO, ^ro) parameter space, 
since Oao can both be negative and positive, so it always 
exists for given flhpo and flro- 

The expression under the root is the discriminant A2 
of the polynomial in equation (|III.16p , when considered 
as quadratic in a;^, i.e., 

A2 - nlpo + r!,o(J^rO - 2) - 2nhpQ{nrQ + 1) + 1. (III.18) 

All signs of A2 are possible for positive values of the 
parameters, and we consider them consecutively. 



A2 < 



In this case the roots Xq^. are complex and different, so 
there are four distinct complex roots of the main polyno- 
mial ljIII.14|l . Accordingly, this is the case B.l.a.l of the 
general classification from Appendix B. The evolution of 
the model is very interesting from the current observa- 
tional point of view, since it starts from Big-Bang and 
terminates at Big- Rip (compare Fig. 



2. A2 > 

Here the roots Xq^ are real, so we analyse their signs, 
which, in turn, depend on the sign of Oao = 1 ~ i^hpo ~ 
flro- It is worth noticing that for both these cases we 
always have 



I^^Aol > 



(III.19) 



For Q\o > we have Vlro < 1 — ^Ihpo, and the roots 
a^oi are negative thanks to the inequality ljIII.19|) . and 
this is the case B.l.a.l of Appendix B again. 

For r^Ao < 0, we have firo > 1 — i^hpo, and now both 
roots Xq_|_ are positive, so the main polynomial has four 
distinct real roots. This corresponds to the general case 
B.l.a.3. As there are two positive and two negative roots, 
a possible oscillating solution is singular at y = (BB- 
BR type). There is also a bouncing Big- Rip solution only 
(cf. Fig. USI). 



3. A2 







The conditions for the parameters are 

n^pO = (1 ± V^of ^ ^rO = (1 ± \A4^)', (111.20) 

which imply that 



- —'^\/^hpo{\/^hpO ± 1) ■ 

(III.21) 

We thus obtain two double roots, given by 



xo± 



AO _ \/ ^hpO ± 1 
\/^hpO 



2n 



hpO 



(III.22) 



If f^Ao < we have two double real roots - one pair 
negative, one positive. That is, the general case B.l.c.l 
so that there are two unstable static (US) solutions. Be- 
sides, there is an interesting asymptotic solution which is 
monotonic from one static solution to the other (MUS- 
MUS type) , and the two bouncing solutions - one of them 
asymptotes from a static towards Big-Bang and another 
from a static towards Big- Rip (cf. Fig. I18|l . 

If f2Ao > we have two double complex roots, which is 
the general case B.l.c.2 - the solution is of BB-BR type 
(cf. Fig. UH). 

Finally, for f^Ao = we obtain a quadruple root a; = 0, 
and the general solution is that of B.l.e. However, since 
an unstable static universe is at y = 0, then only one of 
the asymptotic solutions is physical and is again a BB- 
BR type (cf. Fig. 1221). 



C. Superphantom and radiation models 

Allowing superphantom (7 = —2/3) and radiation 
(7 = 4/3) only in the universe we get from ljIII.l|) and 
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(llIL2l (cf. Ref. gg) 



dy 



dr 



^spoy* + ^roy 



The equation (|111.23|l solves parametrically by 
ao 



a{ri) 



t{r,) = 



spO 



spO 



Ho 



dr] 



(III.23) 



(III.24) 



(III.25) 



The solution is a BB-BR type. The Big-Bang sigularity 
at a — > appears for 



Vbb = T]o 



n 



spO 



(III.26) 



while the Big-Rip singularity at a — > oo appears for 

rj-rjo^oo . (III.27) 

The energy density of radiation and superphantom equal- 
ity time is given by 



Veq 



2m 



spO 



2n 



+ a 



spQ 



tO 



(III.28) 



This is the best exact solution to investigate the proper- 
ties of phantom model evolution from Big-Bang to Big- 
Rip, since it integrates in elementary functions. As it was 
shown in Ref. j45| this is not the case for dust (7 = 0) 
and phantom (7 = —2/3) only model. 



IV. OBSERVATIONAL QUANTITIES FOR 
GENERAL PHANTOM MODELS 

A. 0(2^) luminosity distance Dl{z) formula 

We will follow standard derivation of the luminosity 
distance |54ll. b ut include higher order characteristics 
such as jerk |53| and "kerk" 111.23^ as state-finders. The 
physical distance D which is travelled by a light ray emit- 
ted at the time ii, and received at the time to is given 
by 



D 



to 



dt — c{to — ti) . 



(IV.l) 



The scale factor a{t) and its inverse 1 / a(t) at any moment 
of time t can be obtained as series expansion around to 
as (a(<o) = ao) 

a{t) = ao |l + Ho{t - to) - ^qoH^{t - tof (IV.2) 
+ lijoHUt - tof - l^koH^it - tor + 0[it - tof] 



and 



^ = i^{l + i/o(t„-t)+i.„^(f + l)(to-t)^ 

(IV.3) 



-^0 ( 90 + f + 1 ) {to - 



+0[{to'tf]} , 



Putting t = ti in (|IV.3|I . and using l|II.18|l and l|IV.l(l . we 

get a series z = z{D) as follows 



l + z= 1 



HoD 



1 ^^.(HoD 
-(go + 2)^ 



1 r . M ( HoD 
+ g bo + 6(go + 1)] f-^ 



'7^ [ko + 8jo + 6go(go + 6) + 24] 
24 V c 



-0 



HoD 



(IV.4) 



This series can be inverted to get D = D{z), i.e., 



m = — o - ( 1 



90 \ , /, , ,95 Jo \ 2 
-jz+ll + qo + --j]z 



, jo /5 \ ko 53 3 ■ 

+0{z^)} . (IV.5) 

The luminosity distance Dl is defined as the distance to 
an object whose energy flux falls off in the way, as if we 
dealt with a Euclidean space j54] . i.e., 



F = 



L 



(IV.6) 



where L is the total brightness of the source. In the 
Friedmann universe it reads as 



D, 



(1 + 2;)aoro = (1 + z)aoSK{x)i 



(IV. 7) 



where tq is the radial distance from a source to an ob- 
server, and 

^sin(Vi?x),if >0 

Sk{x)^{ X,K = (iv.8) 
/^sinh(/p^x),^<0 



From the null geodesic equation in the Friedmann uni- 
verse we have 



cdt 



dr 



Vl - Kr^ 



= Xiro) = Sj^^iro) , 



(IV.9) 



8 



so that we can conclude that 



ro = Sk{x) = S_ 



K 



cdt 



(IV. 10) 



/*, a{t)J ' 

which can be expanded in series for small distances as 
/ / cdt \ K I f" cdt \ ^ I cdt ^ ^ 



cdt 



*" cdt\ K 

(iv.ii) 

Now, the problem reduces to a calculation of an integral 
(|IV.9p in terms of the distance D. For this sake, we use 
series expansion (|IV.3|I . i.e., 

cdt __ 

= - / " dJi + i/o(to - + al (? + l) (io - tf 



+O[(io-t)']} , 
which integrates to give 



cdt _ D \ (HqD 

ti ao I 2 V c 



1/90 



3 V 2 



5 V 3 4 



HoD 



2* 24 



01 



Notice that 



1 



'42 



3 /'i?o£' 
2 



c 
c 



(IV. 12) 



(IV. 13) 



(IV. 14) 




1 

12 



KO 



(IV. 15) 

In view of the definition ljIV.7|) . it is useful to express tq 
as a function of redshift z, using formula HIV.5p . i.e.. 



' 6 



6 



5 ko 9 ^33/32 

Y^9oJo-^ + ^Jo- I 290 



90 + 77 



KO 



z^ + Oiz^) 



(IV.16) 



Finally, from ljIV.7|) and ljIV.16p . we have the fourth- 
order in redshift z formula for the luminosity distance 

D,{z) = ^x 
no 

1 + ^(1 - qo)z + \ [90(390 + 1) - (jo + 1) - ^Ko] 



1 

24 



[5jo(29o + 1) - fco - 159o(9o + 1) + 2(1 - qo) 



+2nKo{iqo + l)]z^ + 0{z^)] . (IV.17) 

The formula (IIV.17|I agrees with the third-order in z for- 
mula ( 37) of {52] and with the fourth-order in z formula 
(7) of 1^. It also agrees with the formula (39) of |3| for 

^stO = ^gO = ^phO — ^spO = ^bpO = ^hpO- 

This series expansion for the luminosity distance 
can alternatively be obtained by using the method of 
Kristian and Sachs 2, 58^, 5^ . 

Notice that one can reduce nicely the term of the third- 
order in redshift z by using the formulas ljll.36|) and 
(|II.32|) to get (we dropped the fourth-order term which 
can also be reduced using these expressions together with 
(llL37ll ^ 

Dl 



cz { \, 
— U + -(1 
Ho I 2^ 

1 



6 



9o)z+ 



[390(90 + l)-15f^sto-10f^ 



90 



(IV.18) 



phO 



-6rif,po ~ lOilftpo] z^ 



3^^sp0 
0(z3)} 



Using ^IV.im . <IV.13II . and IITV.14II we get 



B. Luminosity distance - some special cases 



roiP) 



aoHo 



+ i(2 + ,o-f^A'o) 



(90 + ^ + 1 - ^KO 



1 fHoD 
2 



HnD 



Let us now start with the exact luminosity distance 
formula in the two simple cases of only standard matter 
(|II.14I) or phantom 111.151) , respectively, which read as 
(compare 0|) 

Dr. = ,^^^ [(l + z)-(l. 



37 - 2 i/o 



(7 ^ 2/3) , 
(IV. 19) 



9 



and 
Dl 



3 I 7 I +2 i^o 



3|-y|+4 

(l + z)-(l + z)^ 



(IV.20) 



reads as 



Dr ^ —n 



2 — ri^'Q — 2fi 



rO 



spO 



In the special case of Section III. A, if we apply the fact 
that the comoving distance x is related to the cosmolog- 
ical time according to ljIV.9p . and change the variable 
X — (a/ao)^, we obtain 



2Hoao 



2/1 \ 2 

da; 



-) =xW{x), 



(IV.21) 



which can again be transformed into the Weierstrass 
equation, as W is of degree three only. Consequently, 
we can obtain the following 



Ho 



~ 2nK'o{2 - 617rt + i^spo)] + C'(z3)}(JV.26) 

and it is in agreement with our general formula ljIV.18|) 
after rearranging the terms by using (|II. 32^ - 111. 36|) . 



C. 0{z'^) redshift-magnitude relation m{z) formula 

As for the apparent magnitude m of the source we have 



2i7oao , 
P( X) = Px 



(IV.22) 



m = 5 logjQ Dl+ M , 



(IV.27) 



where Af is the absolute magnitude. We will expand 
12z^{2 + zY ' ' ' ' the apparent magnitude into series using the luminosity 

(1 + zf(2 + 2z + z^) + f^A (2 + 2z + z^) + 2n. 1 expansion \lV.n\. Following the Refs. '5^,'5l, 

° '■''^^ in order to adopt to the standard astronomical notation 

for the magnitude we first write (|IV.27|I as 



[nK'oiz^ + /^z? + lOz^ + 12z + 6) 



+ 6i/a-o(l + zf + f7if'o(l + zf + nAo(l + zY + Ospo \ , 



where we have used the definition of redshift (|II.18p . The 
function p is not the same as in an appropriate cosmo- 
logical solution. The invariants here are 



m - M = - logio 



(IV.28) 



92 



— f2^/Q rilAnf^ 

1 



AO"rOj 



Using lIV.lTp . one easily gets 

m~ M = 5 logio(c2:) - 5 log^g -f^o 



y logio{l + (1 - 1a)z 



Substituting HIV.23|) into (|IVJI), we get 

Dl 



15 2 qo 



Hq a 



[jo(10go + 3) 
+0(z4)} . 



KO 



(IV.29) 



3go(5go + 3) + 6qo^Ko] z^ 



where the new subscript of S indicates that instead of K 
we need to use the related density fix- 

Since we know that a; = is necessarily a root of the 
right hand side of equation ljIV.2ip . we can easily obtain 
a special case with a double root, if we put Qro = 0. The 
Weierstrass function will then reduce to a trigonometric 
one, and the luminosity distance will become 



This can further be expanded using the properties of log- 
arithms as 



5 51n5r(z) 



-logioe)lng(z) , (IV.30) 



Dj 



c(l 



Hq 



1 



y^—flK'O 



arccot 



4Fv 



2 

Hk'o 3 

(iv'24) 

If we assume that flcso — 0, or in other words Qk'o = 
— Oif , the above formula is further simplified to 



where g{z) is a given function of redshift z, to get 

'5. 



m- M = 5 logio (cz) - 5 logio Hq - log 



=10 ' 



(1 - qa)z + 



1 



90 
2 



7 

- jo — ^ - ^KO 



Dr 



c(l + z) 



^ [2jo {8qo + 5) - 2fco - qo {Iql + Hqo + 23) + 25 



HoxI'^Px + k^K 



(IV.25) 



+mKo{2qo + l)]z^ + 0{z^)] 



(IV.31) 



For practical application, it might also be useful to con- 
sider the series expansion of the formula ljIV.23p . which 



From (|IV.31|) it is clear that the jerk appears in the sec- 
ond order of the expansion and the "kerk" appears in 
the third order of this expansion. Also, if 0{z^) term 



10 



is drooped, this formula agrees with the formula (41) of 
Ref. P, i.e., 



m = M - 5 logio Ho + 5 log^o (cz) 

1 



+ (2.5 logioe) I (1-90)2 + 
2a sQ 2Aq 



(IV.32) 

(3(7o + l)(go-l)- 



2^ + 0(z^) , 



in which all but radiation, dust, A-term, and cosmic 
string terms are neglected. This can easily be checked 
after the application of the relations ljll.321 (|II.34p . and 
(jll.36p . The most general, second order 0{z'^) relation, 
for all the types of matter admitted in this paper reads 
as 



TO = M - 5 log^o Ho + 5 logio (cz) 

(2.51ogioe)|(l-go)2+ ^ 
10 



(IV.33) 



1 



-n 



gO 



2a 



rO 



-a 



phO 



spQ 



2n 



bpO 



2 V 2 



-n 



sO 



hpO 



0(z3 



D. Angular Diameter 

This quantity is given by: 

Dl ' 



(IV.34) 



where d is the linear size of a given object. Using the 
formula (|IV.9|) we can consider the minimum of 0(z), 
which is a solution of the equation 



Sk{x) 

s'Ax) 



(IV.35) 



In general this is a transcendental relation, which is im- 
possible to solve in a suitably closed form. However, in a 
special case of J7,o = ^sco = 0, this equation admits an 
algebraic solution, as it simplifies to: 



l2P^-2nK + zfW 



(IV.36) 



Finally, the angular size of a galaxy in a one- 
component standard (|II.14|I , and phantom (|II.15p matter 
models are (d-linear size of this galaxy) 



and 



A(9i 



ffo(37-2) 



d{l + zf 



2c 



(l + z)-(l + z)- 



(7 ^ 2/3) , 
(IV.37) 



ffo(3 I 7 I +2) 



2c 



(l + z)-(l + z)^' 



(7 ^ 2/3) 
(IV.38) 



In the former case, there is a minimum at 

2 

2 \ 2-3^ 



37 



1, (7y^2/3); 



(IV.39) 



Zrmn = 6 - 1 , (7 = 2/3) ; Z„„ QO (7 = 0), 

while for the latter (phantom) the minimum does not 
appear at all. 



E. Source Counts 

We can also calculate the number of sources with red- 
shifts from the interval z, z-t-dz and the density of sources 
n(z) 



N = A^n{z)alSl[x{z)]'^, 



can be rewritten using (|IV.21|I as 



N = 47rn(z) 



Ho J {l + z)^^W[{l + zy 



(IV.40) 



(IV.41) 



or, assuming constant n(z) = n — const., it can be ex- 
panded into 



N = 4™ j |z2-2(r!A-'o+21^.rO-f^spo)2'+0(z4) 

(IV.42) 

The O(z^) formula for radiation, matter, strings and A- 
term models was given in as 



/ 2^rO + ff^mO + ^^sO — 90 — 1 



(IV.43) 



1 



X |z^ - 2(90 + l)z-' + — [{qo + l)(379o + 31) 

- 48^2,0 -2m™o-4r!,o] 2^ + 0(z^)}. 

However, both angular diameter and source counts 
tests, in view of supernovae data are not so precise 
as redshift-magnitude relation test, which is the reason 
we do not investigate them in so much detailed way, as 
we do for the redshift-magnitude. 



V. CONCLUSION 

We have studied more general phantom {p < —g) cos- 
mological models (cf. 113), which include all the 

types of both phantom and standard matter from the 
barotropic index range —7/3 < w < 1. 

We have found that there are various interesting possi- 
bilities of the evolution, depending on the matter content 
in the universe. Since phantom cosmologies allow both 
standard Big-Bang (a — > 0, if i — > 0) and phantom-driven 
Big Rip (a — s- oo, if i — > to) singularities, then the set of 



11 



possible types of evolution enlarges. From the observa- 
tional point of view the most interesting is the Big-Bang 
to Big-Rip (BB-BR) type of evolution which is strongly 
supported by the current supernovae data I3| . However, a 
lot of other interesting theoretical options related to this 
are also possible. For example, there exists an unstable 
static universe (US) which, if perturbed, may go into one 
of the two monotonic universes - one of them is mono- 
tonic towards a Big-Bang (MBB type) and another one is 
monotonic towards a Big- Rip (MBR type) . More hybrid 
solution is that there exist two unstable static universes, 
if perturbed, except for MBB and MBR options, there is 
a possibility to have a solution which is monotonic from 
one US solution to the other (MUS-MUS type). In fact, 
one or the two monotonic solutions may have a bounce 
- it is of monotonic non-singular (MNS) type. Stable 
static solutions (SS) are also possible. The standard Big- 
Bang to Big-Bang (BB-BB) type of solutions in phantom 
cosmologies may also be replaced by Big-Rip to Big-Rip 
(BR-BR) type. Finally, there exist a large class of non- 
singular oscillating (ONS) type of solutions - in some 
cases oscillations are allowed for the two different ranges 
of the values of the scale factor. 



We have also studied observational characteristics of 
phantom cosmologies. We enlarged the set of the dynam- 
ical parameters of cosmological models to include (except 
for the Hubble constant and the deceleration parame- 
ter) the jerk and the "kerk". The "kerk" is defined by 
the fourth order derivative of the scale factor. This al- 
lowed us to write down the luminosity distance relation 
Dl{z) and the redshift-magnitude relation m[z) up to 
the fourth order term in the series expansion of redshift 
z. We conclude that jerk appears in the third order of the 
expansion while the "kerk" appears in the fourth order 
of the expansion, as expected. Both jerk and "kerk" can 
be useful as state-finders - the parameters which can be 
helpful in determination of the equation of state of the 
cosmic fluid [s^, . Of course these considerations for 
state-finders are valid, provided we assume a barotropic 
(or analytic) type of the equation of state and do not 
apply in the case of sudden future singularities which do 
not tight p and g [Tllli^ . 



Finally, we have discussed other observational tests for 
phantom models such as angular diameter test and source 
counts. It is interesting that for one fluid phantom mod- 
els the angular diameter minimum does not appear at all 
- this is in clear contrast to a standard matter models 
where the minimum can be an important characteristic 
of these models. 



As a further step towards the more negative values 
of pressure in cosmology (despite some objections [lOjV 
phantom cosmologies are a viable completion of standard 
cosmologies which may solve cosmological puzzles. 
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APPENDIX A: GENERAL CLASSIFICATION OF 

THE MODELS WHICH LEAD TO A CUBIC 
POLYNOMIAL IN THE CANONICAL EQUATION 

We take into account the canonical equation 

— j =azx^ + a2X^ +aix + ao^W{x) (A.l) 



together with the constraint 



03 + 02 -I- ai -I- flo = 1, 



(A.2) 



for all non-static solutions. The reason for this is the 
normalization of an independent variable, so that for a 
particular value of x its derivative is equal to a fixed 
constant. Here we selected W{\) = 1. This allows eas- 
ier classification, as it directly corresponds to a shape 
of the polynomial W{x). Since the parts of the curve 
{{W{x),x) : W{x) > 0} will in general be disconnected, 
a physically significant solution will be the one for which 
X passes through 1. 

In the case of a static solution, W{x) = at all times, 
and such a normalization is impossible. 

The behaviour of the solutions depends primarily on 
the roots of the equation W{x) ~ 0. As we are only 
interested in real and positive solutions, the position of 
the roots, determine the intervals where W{x) is positive, 
i.e., where physically significant evolution can take place. 
To this end, we introduce the following notation: 

Bi : W{x) ^ a3{x ~ ei){x - e2){x - es), 
A := (61-62)^(62-63)^(63-61)^ 



alal al afao 020100 „„o§ 
— 4 ^^"3 ~ ^ — T' + 3 

Ct^ ^3 ^3 ^3 ^"3 



1 2 1 

92 := ^2^2 ~ 4«30i, 



33 := 



48 



030201 



1 3 1 2 



216 ^ 16 



O3O0, 



where are the roots of W{x), A is the discriminant, 
and (72, 33 are invariants [H^ . For static solutions, the 
appropriate quantities are those of equation ljA.l|l before 
it is normalized. 

The invariants (72 and 173 appear when equation IjA.ip 
is transformed by 



4 02 
x = —V - — 
03 Sas 



(A.3) 



It is then simplified to 



di^\' .3 



— J = iv-" ~ g2V~ gs^Viv), (A.4) 
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and can immediately be solved by means of the Weier- 
strass elliptic p function, i.e., 

V ^ p{t - To;g2,g3)- (A.5) 
Now, we consider the following cases: 



as > 



a. A > 



Necessarily, there exist three distinctive, real roots of 
the polynomial W{x) in this case. Arranging them so 
that ei > 62 > 63, the condition (|A.2p implies 63 < 
1 < 62 or 61 < 1. Accordingly, there are two regions 
of admissible values of x: 63 < a;_ < 62 and 61 < a;+, 
respectively. An appropriate general solution is 

_ P'(t - To) + 0321 [p(t - To) - j^(3a3 + 02)] + \a3 



2[p(t - To) - ^(3a3 + as)] 



with 



(A.6) 



0321 



-(3a3 + 2a2 + oi) , 



or, for particular cases 



<13 



63 



63 



p{t - to + c^3) - II 

_ 3a3(e3-ei)(e3-e2) 
12p(r-ro)-(2e3-ei-e2) 

24p(r-ro)-H"'(e3) ' 

-p(t-to)- f^. 

as " ^ 303 



(A.7) 
(A.8) 



Here To is a real constant, and 103 is a purely imaginary 
half-period of p, given by 



UJ3 



dt 



93 



The solution x_ is a non-singular solution, if 63 > 0, 
while x+ is singular, when t — To = 2nuji, n g Z (Fig^. 

Both of the solutions (|A.7|I and l)A.8(l are periodic with 
the real period 



dt 



^4^' - 92t 



93 



b. A = 

Here, there is a double or a triple root. In order to 
distinguish between these two possibilities, we use the 
equation ljA.4|l . The roots of the polynomial V{v): ei, 
62, 63, must satisfy: e 1 -f 62 + 63 — 0. Therefore a triple 
root of W{x) — corresponds to a triple root w = of 
V{v), with 32 = 53 = 0. 



FIG. 1: Case A.l.a - three different roots. There is an oscil- 
lating non-singular (ONS) solution X- given by <A.7ll and a 
singular solution (of Big-Rip to Big-Rip (BR-BR) type) x+ 
given by IIA.8II . 



Let us investigate the former case first. Obviously, 
there must exist another real root, but not necessarily 
positive. That is, W{x) = 03(0; — ei){x — 62)^. 

If the double root is a smaller one, then because 03 > 0, 
we obtain a stable, static solution x — ei. If the opposite 
is true, the same static solution is unstable. 

For 62 < 61 (61 < 1) we have the following solution 
with x>ei (Fig. 0) 



Xo 



ei-l-(62 — 6i)tan^ 



ei - 62) 



(A.9) 



which is periodic with a half-period 

TT 

UJl 



y/a3{ei - 62) 




FIG. 2: Case A.l.b - the smaller root 62 is double. There is a 
stable static solution (SS) x = ei and a solution xo (BR-BR 
type) given by (IA.9t . 

In fact, this is valid for both 62 < 61 and 62 > 61. 
However, when 62 > 61, the solution is no longer periodic, 
and, upon eliminating the imaginary part of To, it is more 
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convenient to write it as (Fig. O 



X- — ei + (e2 — ei) tanh^ 
when ei < 1 < 62 

x+ = ei + (e2 — ei) coth^ 
when 62 < 1. 



^03(62 - 61) 



^03(62 - 61) 



(A.IO) 
(t-t-q) , 
(A.ll) 




FIG. 4: Case A.l.b - one triple root. There is an unstable 
static solution (US) x — e2 > and a monotonic solution 
(MBR) xo given by 11021 . 



FIG. 3: Case A.l.b - the greater root 62 is double. There is an 
unstable static solution (US) x — 62 and the two monotonic 
solutions x± given by lA.lOII and lA.lH . One of them is 
monotonic non-singular (MNS) and another monotonic to a 
Big- Rip singularity (MBR). 

In terms of the coefficients of the polynomial W{x), it 
is possible to discern between these two possibilities, as 
clearly seen from the sign of W{xi). The value at the 
point of inflection is positive when 62 > 61, and nega- 
tive otherwise. This relation can be precisely written as: 
1633 = -alW{xi) = ^ai(ei - ea)^ 

If there is only one triple root 62, there is an unstable 
static solution x = €2, and an asymptotic solution (Fig. 
EI) 



Xq ^ €2 + 



(A.12) 



That is, a "small perturbation" causes the monotonic 
evolution. Of course, if 62 < 0, only the latter solution 
remains valid and it is of BB-BR type. 



c. A < 

Only one real root exists in this case, and because of 
(IA.2p . it must be smaller than unity. Then, there is only 
one type of solution 



4 ^ 02 

Xo = —P(T - To) - - — , 
03 3a3 



(A.13) 



which is singular when the argument r — tq 
4r^7^e(wl), n G Z (Fig. EJ. 




FIG. 5: Case A.l.c - one real root exists only. The solution 
Xo is given by IIA.13t . It is a BR-BR type of solution with no 
Big-Bang (BB) singularity. 



2. a:i<0 
a. A > 

This case is, in fact, the same as 03 > 0, with an 
interchange of possible regions of the evolution. This 
gives rise to a periodic solution oscillating between the 
two greater roots 62 < x < 61 (Fig. 



X+ = 61 



303(61 - 62)(6i - 63) 

12p(r - To) - (261 - 62 - 63) ' 



(A.14) 



where the roots have been arranged as before. 

Another solution, which in the previous case A.l was 
singular, now is also bounded < x < 63, and is given 
by (Fig. El 



4 s a2 

— p(r-To) - — . 
03 3a3 



(A.15) 
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static (US) solution x = 62- The other two describe the 
motion between < x < 62 (asymptotic to the previous 
one) and 62 < x < ei. They are given by (Fig. ISJ 



X- = ei + (e2 — ei) coth^ 

when 1 < 62 
x+ = ei + (e2 — ei) tanh^ 

when 62 < 1 < ei. 



v/|a3|(ei - 62) 



VlasKei - 62) 



(t'-t-o) , 
(A.17) 

(t-to) , 
(A.18) 



FIG. 6: Case A. 2. a - three different roots. Tliere is an oscil- 
lating non-singular solution (ONS) x+ given by 1A.14II . and a 
singular solution (of Big-Bang to Big-Crunch (BB-BC) type) 
X- given by 1A.15II . 



b. A = 

As in the case A.l.b of the 03 > case, we can either 
have a double or a triple root. Also, we can distinguish 
between the two double root subcases, using the same 
formulas as before. 

If the root is only double we have a situation similar 
to what we had before. It can be a greater root, in which 
case the other root must satisfy ei > 1, for only then we 
could have W{1) — 1. We can only have a stable, static 
solution X — 62, and another one for < x < ei, given 



by (Fig.CI): 



xo = ei + (ei — 62) tan^ 



V|a3|(e2 - ei) 



(r-ro) 



(A.16) 

(Again, this is valid for all double-root solutions, but tq 
might be complex.) 



FIG. 8: Case A.2.b - the smaller root 62 is double. There is 
an unstable static (US) solution x = ei and the two mono- 
tonic solutions x± given by ljA.17^ and ljA.18|l . The former 
is of monotonia non-singular (MNS) type and the latter is 
monotonic Big-Bang type (MBB). 

When 62 is a triple root, there is an unstable, static 
solution a; = 62, and a monotonic solution in the form 

(Fig. m 



Xq = e2 + 



asir - Toy 



(A.19) 



In order for condition (|A.2|I to be satisfied, we must 
have 62 > 1, as W{x) is positive only for x < 62. 




FIG. 7: Case A.2.b - the greater root 62 is double. There is a 
stable static solution (SS) x = ei and a BB-BC type singular 
solution xo given by IIA.16t . 

On the other hand, when 62 < ei, which requires ei > 
1, there are three possible solutions. One is an unstable 



c. A < 

As in the case l.c, there is only one real root 61, only 
this time it limits the values of a; to < a; < 61. The 
solution is singular of Big-Bang to Big-Crunch (BB-BC) 
type (Fig. ^ 



4 ,02 
Xo = —p[r - To) - — , 
03 6a^ 



(A.20) 



3. as = 

This case was studied in the Appendix A of Ref. ^£ 
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FIG. 9: Case A.2.b - one triple root. There is an unstable 
static (US) solution x = e2 and a monotonic Big-Bang (MBB) 
solution 2:0 given by lUlOl . 




FIG. 10: Case 2.c - one real root only. The solution xq is 
given by and is of BB-BC type. 



APPENDIX B: GENERAL CLASSIFICATION OF 

THE MODELS WHICH LEAD TO A QUARTIC 
POLYNOMIAL IN THE CANONICAL EQUATION 

In this Appendix we consider the canonical equation 

/da;\ ^ 

( — 1 = a^x'^ + a'iX^ + a2x'^ + aix ao — W {x) , (B.l) 

and impose the constraint 

W^(l) = 04 + 03 + a2 + ai +ao = 1 (B.2) 

which holds with the same restrictions as in the previous 
case of Appendix A. In order to make the analysis of the 
roots easier, we transform (|B.2(I to the form 



— W{x—- — )=x + px + qx + r. 
4a4 



(B.3) 



Because of that, it is easier to define a set of auxiliary 
quantities 



where are the four roots of the polynomial W{x), and 
i,j = 1,2,3,4. It is clear, that the original equation 
(|B.lll . and the transformed equation l|B.3(l have exactly 
the same values of S, and general properties of the roots 
remain the same. Now let be the set of all possible 
values of the symmetric multi-index ij. We can further 
define the quantities: 



^6 = 



CT5 



0-4 = 



(^3 



0^2 



0-1 



E n 



Kl,K2eJ',Kl<K2 

E 



^K2 ^K-i 1 



E 

^1 JT'AI <'^2 



(B.4) 



As these quantities are symmetric in the roots, they can 
be re-expressed by using the coefficients or, equiva- 
lently, p, q and r as 



CT6 = -4pV - 27g'* -f 16/ 



lAApq^r ~ 128pV^ 



CT5 = -4p'' - 18p^g^ 
0-4 



32p'^r - 2lQq^r + 192pr^ 



17/ + A8pq^ 



2-ip^r 



112r^ 



(73 = -28p^ - 2Qq^ - 16pr, 



256r^ 



(B.5) 



(72 = 22/ 

(71 = -8p 



8r, 



Note that erg can immediately be identified as the origi- 
nal equation's discriminant. Together with other cr's, it 
determines the behaviour of the roots in the following 
way: 

1. ^ 

There are four distinct roots, with the following 
subcases: 

(76 < - two complex, conjugate roots, and two 
real roots, 

(76 > - four real roots, or four complex roots, 
CO jugate in pairs. 

2. (76 = 0, (75 ^ 

There is one double, real root, and: 
(75 > - two real roots, 
(75 < - two complex conjugate roots. 

3. (76 = = (75 = 0, (74 7^ 

There are two double roots, which for: 
(7i > are real, 

(7i < are complex, conjugate in pairs. 

4. (76 = (75 = (74 = 0, (73 ^ 

There is a triple real root. The fourth root is also 
real. 
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5. (Te = ""s = "'4 = 0, (73 = <J4> cr2 = (Ti = 
There is a quadruple root - it is necessarily real. 

In order to distinguish between the real and complex 
roots when > 0, we choose to employ the resolvent 
polynomial of W{x), which appears in the process of fac- 
torizing a quartic polynomial into two quadratic ones 

+ px^ + qx + r ~ (x^ + kx + m){x^ — kx + n). 

It is essentially the equation determining the value of fc, 
as for w ~ k'^, the resolvent equation is 

+ 2pw'^ + {p^ - 4r)w - = 0. 

Now, in the case of the four real roots, such decomposi- 
tion is possible in three different ways, as there are three 
ways of pairing the roots. Hence, there must be three, 
distinct, positive solutions for fc^. If, however, the roots 
are all complex, in order for the coefficients fc, m, n to be 
real, each root must be paired with its conjugate. Thus, 
there is only one positive solution for /c^, and the two 
negative ones. It is easy to distinguish between these 
possibilities, looking at the positions of the extrema of 
the resolvent polynomial. In the first case, there must 
be the two positive ones, and in the second, if any exist, 
at least one must be negative. Thus, the problem may 
be reduced to an investigation of the sign of the smallest 
extremum: 

= i(-2p- v'p^-M2r). (B.6) 

If p^ + 12r is negative, and there are no extrema, only 
one decomposition is possible. Note that this quantity 
must be non-negative in all the multiple root cases, as 
degeneracy of the roots of W implies multiple roots of 
the resolvent. 

In general, when there are no multiple roots, the main 
equation ljB.l|l . can be solved by using the substitution: 

1 _ Av W"{e,) 



x-e, W'{ei) 6VF'(e,)' 
which transforms it into the Weierstrass equation 

with the following invariants: 



(B.7) 



1 



92 



1 

12' 



93 = ^040200 

6 



16 



2 



48 
1 



a3a2ai 



1 2 1 
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2, 



-pr). 



216 
(B.9) 



16" 6 

This, however, requires the knowledge of the roots, and 
usually yields too cumbersome formulas, sometimes with 
explicitly imaginary coefficients. The following solutions 
were possible to be simplified as a result of the division 
into the special subcases. Still, the main function p(t), 
and its half-periods, are constructed using the invariants 
(iRflll. 



a4 > 



a. (T6 7^ - simple roots 



Choosing an integration constant by imposing that 
x(0) = 1, we obtain a general solution of the form 



= 1 + 



(B.IO) 



p'(T-To) + i(4a4 + a3) 



2[p(r - To) - 3^(604 -I- 3a3 + 02)]^ - ^04 
i(4a4 + 3a3 -I- 2a2 + al)[p(r - tq) - 5^(604 + 803 + 02)] 
2[p(t - To) - j^(6a4 -f 803 + 02)]^ - ia4 

which is valid for all subcases. However, the behaviour 
of this solution varies greatly depending on the proper- 
ties of the roots. The details are given in the following 
subsections, and are clearly drawn in the figures. 

a.l. Four complex roots 

The polynomial W is always positive here, and from a 
simple, geometric investigation it is clear that the condi- 
tion (|B.2p can be satisfied. Only one solution is present, 
with a possible time reversal t —t applicable in all 
cases. It reaches both Big-Bang and Big-Rip in finite 
times (Fig. [TT|) . 



FIG. 11: Case B.l.a.l - four complex roots. The solution 
xo is given by the formula IIB.IOII - the evolution starts from 
Big-Bang and terminates at Big- Rip (BR-BR type). 



a. 2. Two complex and two real roots 

With real roots, the possible domain of x is separated 
into two regions: a;_ < ei < 62 < x^. However, both 

If 



solutions are given by the same generic formula IjB.lOl 
we are interested only in physical solutions, then only one 
"branch" is valid in each case, depending on the position 
of the roots. Namely, if 1 < ei, we take a;_, and if 
62 < 1, we take x+ (Fig. I12|l . The solution x- starts at 
Big-Bang reaches the maximum and terminates at Big- 
Crunch. The solution x+ starts at Big-Rip reaches the 
minimum and terminates at Big-Rip again. 
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FIG. 12: Case B.l.a.2 - two complex and two real roots. 
There are two solutions x± both given by general formula 
HB.IOII . The solution X- is of BB-BC type while the the so- 
lution x+ is of BR-BR type. 



a. 3. Four real roots 

The situation is similar to the case a. 2 and with two 
more roots, there is a third admissible region only. We 
now have: a;„ < ei < 62 < xq < 63 < 64 < x+, depend- 
ing on which of these intervals the line x = I belongs to. 
The solutions x_ (BB to BC) and x+ (BR only), are de- 
scribed by equation (|B.10I) . and xq is obtained, by adding 
a purely imaginary half-period W3 to the argument. This 
last solution is oscillating with a period 2uji, and avoids 
any singularity when < 62 (Fig. I13|l . 



In this case, the double root is real, and the remaining 
two roots are complex. Using equation (|B.3() . we can 
easily obtain a double root 



ei 



-3q 



4:p + 2y/jP+T2r 



0-3 

404' 



(B.ll) 



which applied gives an appropriate solution 



x± 



wi" 



3W{'wi'^^ - W{"'^ sinh 



hWi'ir to) 
(B.12) 

where the index indicates that the value W{x) should be 
taken at x = ei. As in l.a.2, this formula incorporates 
both x+ > ei and x_ < ei. Which of these should be 
used depends on whether ei < 1 or 1 < ei, respectively. 
Also, there is an unstable static solution xq = e\ and the 
solutions x± which approach it from either Big-Bang or 
Big-Rip (which differs from standard case of an Einstein 
Static Universe - Fig. ITU) . 




FIG. 13: Case B.l.a.3 
solutions: (oscillating 
by ifBlnt and ||rT1 



four real roots. There are three 
and x± (BR only, BB to BC) given 
respectively. 



- one double root 



b.l. (J5 < - two real and two complex roots 

With the appearance of multiple roots, the solutions 
simplify significantly. Firstly, it is possible to express 
the roots themselves in a suitably short form. Secondly, 
the elliptic function p reduces to a trigonometric or a 
hyperbolic function. 



FIG. 14: Case B.l. b.l - two real and two complex roots. 
There is an unstable static solution x = ei and the two asymp- 
totic solutions x± given by lB.12t . One of them approaches 
Big-Bang (MBB type) and another approaches Big-Rip (MBR 
type). 



h.2. 



> - four real roots 



Depending on the relative values of the roots, three 
situations are possible. Denoting a double root by ei, 
given by the same formula as before, we always have a 
stable static solution xq = ei. Since we know ei, it is also 
possible to obtain the remaining roots, and check which 
of the following is applicable. If 62 < ei < 63, 



X± = 62 



(B.13) 



(62 - 6i)(62 - 63) 



ei 



62 + (ei - 63) tan 1^04(61 - 62)(63 - ei)(T - tq) 



with the choice of the -f or - "branch" depending on 
whether 1 < 62, or 63 < 1 (Fig. [T3|l . 

The second possibility is that 61 < 62 < 63, where the 
solutions in the regions adjacent to 61, that is, < x < 61 
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(62 - ei)(e2 - 63) 



ei - 62 + (eg - ei)coth i^a4(ei - e2)(ei - e3)(r - tq) 
a;_i,2=e2+ (B.17) 

(62 - 6i)(e2 - 63) 



61 - 62 + (63 - ei) tanh ^^04(61 - e2)(ei - 63)(t - to) 

and the static solution xq — e\ is also unstable (Fig. I17|l . 
Here we have 



^a4(ei - 62)(6i - 63) 



zartanh 



61 - 62 
61 - 63 



FIG. 15: Case B.l.b.2 - the double root is between the other 
two roots. There is a stable static solution x = e\ and the 
two solutions given by IIB.13t (of BB-BC and BR-BR type). 



and 61 < a; < 62, tend to 61 asymptotically, and xq = 61 
becomes unstable. The respective formulas are (Fig. I16f) 



62 



(B.14) 



(62 - 6l)(62 - 63) 



61 - 62 + (63 - 61) tanh 5^^4(62 - 6i)(e3 - ei)(7 
a;-i.2 = 62 + 



^0) 
(B.15) 



(62 - 6i)(62 - 63) 




2 ' — — — — — 



61 - 62 + (63 - 61) coth \\/a/^(e2 - 61) (63 - 61) (r - To) 
For a finite value of r, x reaches infinity and 

2 



Tq—VJ 



1/04(61 - 62) (61 - 63) 



zartanh 



61 - 63 
61 - 62 



FIG. 17: Case B.l.b.2 - the greatest root is double. There is 
an unstable static solution x = e\ and the other solutions 2;+ 
(oscillating) and a;_i,2 (from BB to BC, bounce to BR) given 
by lB.17t and 1B.18II . respectively. 




FIG. 16: Case B.l.b.2 - the smallest root is double. There 
is an unstable static (US) solution x = ei, and the three 
oth er solu tions x+ ( MUS) and x_i,2 (MBB, BR-BR) given 
by IB.lSt and 1B.16II . respectively. 



Lastly, when 62 < 63 < 61, the situation is similar to 
the previous one, only the asymptotic branch of x-\ is 
now the upper one, with 

a;+ = 62 + 



c. (76 = 0-5 = 0, (T4 7^ - two double roots 

c.l. ai > - four real roots 
The formula for the roots is further simplified to 
as / 1 



61,2 



4a4 



(B.18) 



Here —2p is equal to the discriminant of the equation {x— 
ei)(x — 62), and is always positive. As for the solutions, 
we have two static, unstable solutions corresponding to 
= 61.2, and the remaining two are given by (Fig. I18|) 



61 



= 61 



tanh 



coth 



-2^^«4(t - To) 



--pa4(r - To) 



,(B.19) 



(B.20) 



(B.16) 
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FIG. 18: Case B.l.c.l - four real roots, two of them are dou- 
ble. There are two unstable static solutions xo ~ ei,xo = £2 
and the three other solutions x+ and X-1^2 given by HB.IQII 
and IB.2OI 1. respectively. The solution x+ asymptotes from 
one unstable static solution to the other (MUS-MUS type) 
while x-i is of MBR type and X-2 is of MBB type. 



C.2. ai < - four complex roots 

Similarly to l.a.l, the polynomial W is always positive, 
but the solution is simplified to 



, /I ^- 
xo = h \ -p tan 

4a4 V 2 



^pa4{T~To) 



(B.21) 



with p always positive this time (Fig. I19|) . 



There are two subcases, with ei < 62 and ei > 62, re- 
spectively, but the solution for both "branches" in each 
subcase is given by the same formula: 



x± 



0461(61 - 62)^(r - ToY 



462 



04(61 - 62)2(r - ro)2 - 4 

Infinite x is reached for a finite value of t: 

2 



(B.24) 



T — To = tU = — . 

V|a4||ei - 62I 

There is also the unstable, static solution xq — ei (Figs. 




FIG. 20: Case B.l.d - one triple root which is greater than a 
single root. There is an unstable static solution xq = ei and 
the other solutions i± (BB-BC and MBR type) are given by 




FIG. 19: Case B.l.c.2 - four complex roots. The solution xq 
is given by IIB.211 and is of BB-BR type. 



d. (Tg = (T5 = (T4 = 0, (73 7^ - one triple root 



FIG. 21: Case B.l.d - one triple root which is smaller than a 
single root. There is an unstable static solution xo — ei and 
the other two solutions x± (BR to BR, BB only) are given by 



The multiple root now becomes: 



ei = 

and the fourth root is: 
62 = 



4a4 2^'^' 



0-3 

4a4 



(B.22) 



(B.23) 



e. ae — — a4 — ai = - one quadruple root 

This is the simplest case, immediately integrable to: 

1 



X± = 61 



aiir - To) ' 



(B.25) 
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where the root is: 



ei 



0-3 

4a4 ' 



(B.26) 



Again, we have two "branches" separated by the static, 
unstable solution a;o — ei (Fig. 



FIG. 22: Case B.l.e - one quadrupole root. There is an 
unstable static (US) solution and the solutions x± given by 
~725l They are of MBR and MBB type. 



a4 < 




FIG. 23: Case B.2.a.2 - two complex and two real roots. This 
is an oscillating non-singular (ONS) solution xo, given by the 
formula J^Jot. 





a. ae ^ - simple roots 
a.l. Four complex roots 

This case is clearly impossible, as the polynomial W is 
everywhere negative. 

a. 2. Two complex and two real roots 

There is only one possible solution, as W{x) > when 
ei > a: > 62, where ei, 62 are the real roots. The general 
formula IjB.lOfl still holds, but the behaviour of the solu- 
tion is qualitatively different. It has a real period of 2oji, 
defined as before, and avoids any singularity, provided 
that ei > (Fig. [^ . 

a. 3. Four real roots 

There are only two possible solutions when ei < 1 < 62 
or 63 < 1 < 64. They both are also given by the general 
formula (|B.10|I . independently in each case. They are 
oscillating, and moreover, and have the same period 2uji, 
so that, if < ei no singularity is present at all (Fig. 

123. 



b. (T6 = - one double root 



b.l. (75 < - two real and two complex roots 



FIG. 24: Case B.2.a.3 - four real roots. This is a "double" 
oscillating case, which means we have the two oscillating solu- 
tions x± of the same period 2iOi given by the formula ljB.12^ . 
They are non-singular provided ei > 0. 



There is only one possible solution: xq = ei, with the 
root given by (|B.11|I . As W is negative everywhere else, 
this solution is stable. 

b.2. CT5 > - four real roots 

Here, there are two subcases, both with the static so- 
lution Xq = ei. The first one occurs when 62 < ei < £3, 
and there are two other solutions asymptotic to xq (Fig. 



113 



62 



(B.27) 



(e2 - ei)(e2 - 63) 



ei - 62 + (63 - ei) coth ^y^\a4\{ei - 62)(63 - ei)(r - tq) 



62 



(B.28) 



(62 - 6i)(e2 - 63) 



61 - 62 -f (63 - ei)tanh i\/|a4|(ei - e2)(e3 - ei)(r - tq) 

If ei < 62 or 63 < ei, xq becomes stable, and there 
is only one other oscillating solution between 62 and 63 



2 ■ 
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(Fig. EH) 



(62 - ei)(e2 - £3) 



ei - 62 




rf. aa — as = a4 — 0, (73 ^ - one triple root 

(B.29) ggggjj^jg^jjy ^}jjg jg j^g^^ Qj^g case, but depending on the 
position of the triple root, the non-static solution 



+ (61-63) tan 5 V|a4|(6l - 62)(6l - 63)(t - Tq) 



x± 



462 - a46l(6l - 62)^(t -~ Tq)'^ 

4 - 04(61 - 62)2(t - ro)2 



(B.30) 



has either a maximum, for 61 < 62, or a minimum, for 
61 > 62. Also, xq — 61 is an unstable static solution 
(Figs. HTJm . 




FIG. 25: Case B.2.b.2 - one double root which is between the 
two single roots. There is an unstable static (US) solution 
xo = ei and the two asymptotic solutions x± given by IB. 2711 
and fn^Sl. 




FIG. 27: Case B.2.d - one triple root greater than a single 
root. There is an unstable static (US) solution xq — ei and a 
monotonia solution x+ given by l)B.30|l which has a minimum. 



FIG. 26: Case B.2.b.2 - one double root greater than the two 
single roots. There is a stable static (SS) solution xq — ei 
and an oscillating solution x+ (ONS type) given by IIB.29L 



c"6 = o"5 = 0, (T4 7^ - two double roots 



FIG. 28: Case B.2.d - one triple root smaller that a single 
root. There is an unstable static solution xo = ei and a 
monotonic solution x+ given by HB.SOII which has a maximum. 



c.l. ai > - four real roots 

The only possible solutions here, are two static, stable 
ones: Xq = 61^2- The roots are given by formula (|B.18p . e. ere = 0-5 = 0-4 = cri = - one quadruple root 

C.2. CTi < - four complex roots There is only one point where W is not negative, thus 

yielding just the stable, static solution of xq = ei^ where 
W is negative for all x, so this case is impossible. 61 is given by the formula ljB.26|) . 
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3. a4 = 

This case is essentially reduced to a cubic polynomial 
case considered in the Appendix A. 
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